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ABSTRACT 


Oceanographic  ob.arvations  near  the  Chevron  spilling  well  off  the  Mississippi 
Delta  in  March,  1970,  revealed  relative  roles  of  varin"  physical  factors  of  the 
regional  estuarine  system  in  the  behavior  of  oil  slicks.  Surface  stress  from  the 
wind  was  most  Important;  at  speeds  above  1$  mph  the  slick  orientation  was  generally 
determined  by  the  wind  direction.  The  wind  a.l  jo  indirectly  affected  oil  which  was 
sunk  by  dispersant  in  that  wind  wavea  promoted  mixing,  which  in  turn  affected  the 
vertical  stability,  hence  eventually  the  velocity  profile.  Wind  setups  and  setdowns 
were  correlated  with  downward  and  upward  isopycnal  movements,  respectively.  Both 
calculations  and  observations  showed  that  tidal  currents  produced  an  L-shaped  slick 
geometry  when  winds  were  below  about  15  mph.  The  diurnal  rotation  of  the  tidal 
currents  served  to  limit  the  excursion  length  of  oil  frem  the  source,  keeping  it 
short  of  the  nearest  shore.  The  presence  of  iresh  water  from  the  Mississippi  River 
in  the  surface  layer  and  the  consequent  development  of  convergence  lines  often 
formed  a  natural  barrier,  preventing  oil  from  encroaching  upon  the  shore.  Theo¬ 
retical  analysis  using  turbulent  diffusion  theory  disclosed  that  the  area  and  length 
of  a  steady-state  oil  slick  Increased  with  oil  discharge  but  decreased  with  the 
current  speed  and  the  lateral  diffusion  coefficient.  Slick  width  increased  with 
the  rate  of  discharge  and  decreased  with  current  speed  but  was  independent  of  the 
diffusion  coefficient.  The  width/length  ratio  of  an  oil  slick,  which  was  notably 
indipendent  of  current  speed,  was  controlled  only  by  the  diffusion  coefficient  and 
the  oil  discharge.  In  order  to  be  able  to  predict  oil  .slick  behavior,  understanding 
of  the  interaction  between  these  multiple  factors  Is  essential  and,  for  that  end, 
the  use  of  a  multiple-sensor  array  with  automatic  data  transfer  and  processing 
capability  is  recommended. 
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NOTATION 


A  slick  area 

C  concentration;  C(y  -  R)  is  the  concentration  at  the  visible  boundary 

of  the  slick 

F  fetch 

g  acceleration  because  of  gravity;  or  grams 

h  thickness  or  depth 

significant  height  of  waves 

K  diffusion  coefficient,  diffusion  coefficient  of  oil 

K  modified  Bessel  function  of  2nd  Kind  and  Zero  Order 

o 

*  -  -  - . - 

1  scale  length  or  average  eddy  sire 

L  slick  length 

Q  rate  of  emission  (of  oil) 

r  distance  from  the  origin 

R(£)  the  Lagrangian  autocorrelation  function 

t  various  time  parameters 

T  time 

U  mean  horizontal  current  speed 

llw  wind  speed  at  10-meter  height 

v'  turbulent  velocity  in  y  Erection 

2  1/2 

(v'  )  intensity  of  turbulence,  usually  of  oil 


W 

6 

M 

F 

P 


slick  width 

one-half  the  initial  angle  of  dispersion  of  the  slick, 

„  do  .~.l/2 
tan  0  *  —  ■  (v  )  /U 
dx 

millimicrons 
the  lag  time 
density 

the  variance  of  particle  displacements  from  the  source 
in  polar  coordinates,  the  argument 


INTRODUCTION 


When  the  Chevron  production  platform  MP-41C  caught  fire  and  burned  in 
February  1970,  a  situation  was  created  which  at  the  time  appeared  to  be 
potentially  a  prelude  to  a  disaster  of  major  scale.  Efforts  to  extinguish 
the  costly  fire  were  almost  certain  to  lead  to  large  losses  of  crude  oil 
Into  the  sea,  and  much  of  this  seemed  likely  to  wash  ashore  in  various  areas 
of  the  marshes  on  the  mainland. 

Figure  1  indicates  the  location  of  the  accident  (platform  41C) ;  it  is 
north  and  east  of  the  Mississippi  Delta  and  near  the  entrance  of  the  Gulf 
outlet  channel.  Water  depths  in  the  general  vicinity  are  from  20  to  50  feet. 
Before  direct  measurements  were  begun,  a  search  for  literature  giving  indi¬ 
cations  of  circulation  patterns  in  the  area  was  made. 

Scruton  (1956)  measured  24-hour  diurnal  current  profiles  at  several 
stations  in  the  Main  Pass  area.  Two  of  these  were  in  the  passages  into 
Breton  Sound  lying  to  the  east  and  west  of  Breton  Island.  In  these,  ebb 
and  flood  current  direction  were  necessarily  constrained  to  inward  and 
cjutward  flow.  Surface  currents  showed  a  net  outward  flow  at  the  time  of 
measurement,  which  was  in  agreement  with  the  wind  direction  at  the  time. 
Another  station,  away  from  the  passes  and  within  3  miles  of  the  site  of  the 
4 1C  platform,  showed  more  variability  in  flow  direction  but  indicated  even 
higher  residual  outward  (in  this  case  southeastward)  flow  at  the  surface. 
Scruton  felt  that  the  residual  outward  flow  in  this  last  case  was  also 
caused  in  part  by  the  wind  and  by  transfer  of  fresh  water  from  the  river 
out  to  sea  at  the  surface.  He  indicated  that  maximum  currents  at  times  of 
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gure  1.  Location  map  with  stations  occupied  in  study.  A  is  wind  rose  for  long-term  March 
servations.  B  is  wxnd  rose  for  March  1970  from  lightship  data. 


tropic  tides  could  be  1-2  knots  and  that  at  times  of  equatorial  tides  currents 
were  much  lower. 

The  general  hydrography  of  the  area  is  Indicated  on  Figure  1  and  is 
given  in  much  greater  detail  on  chart  1270  of  the  U.S.  Coast  and  Geodetic 
Survey.  An  important  alteration  in  the  area  since  the  time  of  Scruton's 
study  was  the  dredging  of  the  Mississippi  River  Gulf  outlet  channel  to  a 
controlled  depth  maintained  at  36  feet,  with  deposition  of  spoil  along  the 
southwest  side  of  the  channel. 

Orton  (1964)  gave  wind  data  for  subsquare  55  of  Marsden  square  81,  in 
which  the  41C  platform  lies.  For  the  "average"  March,  on  the  basis  of  5,331 
observations,  the  winds  are  indicated  on  wind  rose  A  of  Figure  1.  l.ie  strong 
northwestward  and  westward  components  seemingly  would  favor  drift,  insofar 
as  it  is  wind  controlled,  toward  the  Mississippi  Delta  and  the  marshes  sur¬ 
rounding  Breton  Sound.  The  Weather  Bureau  (1959)  also  has  published  data 
indicating  that  the  frequency  distribution  of  direction  of  wave  approach 
for  winter  and  spring  months  shows  a  predominance  of  waves  coming  from  the 
east  and  southeast  in  the  area  of  the  41C  platform. 

Tidal  predictions  given  in  Coast  and  Geodetic  Survey  (1970)  tables 
for  Breton  Island  for  the  month  of  March  1970  are  shown  in  Figure  2. 

Another  very  important  part  of  the  physical  background  necessary  to 
an  understanding  of  the  situation  is  the  influence  of  waters  discharged 
by  the  Mississippi  River.  Figure  3  gives  the  stage  height  at  New  Orleans 
and  discharge  at  Vicksburg,  Mississippi,  averaged  over  several  years. 

March  is  typically  a  time  of  high  stage  and  discharge,  and  1970  was  no 
exception.  The  discharged  floodwaters  are  highly  turbid  and  cool  (10°C) 
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Figure  2*  Tidal  range  for  Breton  Island  and  wind  speed  and  direction 
observed  at  lightship  during  March  1970. 


and  contrast  sharply  with  the  much  less  turbid  and  generally  warmer  (20eC) 
saline  Gulf  waters. 

Except  for  the  above-mentioned  material,  there  was  little  background 
information  useful  for  judging  what  might  happen  when  oil  began  to  spill  from 
the  Chevron  platform.  For  this  reason  it  was  decided  that,  in  light  of  the 
possibly  dangerous  situation  that  could  develop,  some  basic  oceanographic 
information  should  be  collected  before  the  fire  was  put  out  and  during  the 
spill.  The  Eighth  District  of  the  U.S.  Coast  Guard  readily  agreed  to 
cooperate,  and  measurements  were  begun  on  March  5. 

The  most  important  consideration  was  current,  particularly  surface 
current.  Temperature  and  salinity  were  found  to  be  useful  discriminators 
of  the  riverine  c.\J  marine  water  masses  which  crossed  the  accident  area. 

Wind  records  were  also  maintained  to  establish  the  relation  these  might 
have  to  surface  currents.  Currents  at  depth  were  monitored  inasmuch  as  it 
seemed  possible  that  a  decision  to  use  dispersants  to  emulsify  the  oil 
might  lead  to  spread  of  oil  throughout  the  water  mass.  Tides  were 
measured  directly  since  it  was  not  certain  how  well  the  predicted  curve 
would  agree  with  the  real  situation.  Wave  height  and  wave  period  seemed 
useful  indices  that  might  affect  diffusion  of  the  drifting  oil.  Another 
category  of  observations  was  the  direct  delineation  of  the  orientation  and 
shape  of  the  oil  slicks  in  space  and  through  time. 

In  all,  observations  were  made,  from  four  different  Coast  Guard 
vessels  at  nine  different  anchorages  for  continuous  periods  ranging  from 
5  days  to  6  hours.  Slick  observations  were  made  from  helicopters  on  a 
routine  basis  by  Coast  Guard  personnel,  usually  twice  a  day  in  clear 
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weather.  Overflights  were  made  by  several  different  groups  using  various 
conventional  and  mere  novel  means  of  remote  sensing  to  document  the  extent 
of  oil  slicks.  These  data  have  been  made  available  to  us  by  the  U.S. 
Geological  Survey  and  by  Remote  Sensing  Inc.,  Houston,  Texas.  lhe  National 
Aeronautics  and  Space  Administration  also  made  two  overflights  during  which 
the  C.S.I.  group  obtained  ground  truth  rvations.  This  material  has  not 
yet  been  made  available;  however,  one  of  us  was  able  to  inspect  briefly  the 
imagery  at  the  Manned  Spacecraft  Center  in  Houston.  A  tide  and  wave  gage 
was  suspended  from  Chevron's  MP-41M  rig,  which  is  about  1  mile  from  the 
damaged  structure.  A  bottui.  current  ins Laliat ion  was  also  made  at  an 
additional  site,  but  this  instrument  has  not  yet  been  recovered. 

Complete  records  of  all  temperature,  salinity,  and  current  profiles 
made  from  all  nine  stations  are  presented  in  the  appendix.  Graphical 
records  of  these  and  other  observations  are  presented  throughout  the 
remainder  of  this  report. 

The  41C  platform  caught  fire  at  3  a.m.  on  Tuesday,  February  10,  1970. 
While  it  was  burning  unchecked  most  of  the  escaping  oil,  estimated  by 
Chevron  at  1,000  barrels  per  day,  was  burned  and  did  not  reach  the  water 
surface.  The  work  near  the  platform  necessary  for  quenching  the  fire 
required  pumping  water  onto  the  flame,  however,  and  this  tended  to  wash 
oil  into  the  water,  resulting  in  light  slicks  for  about  a  month  before  the 
fire  was  out.  The  first  attempt  (March  8)  to  extinguish  the  flame  was 
successful,  but  spontaneous  reignition  occurred  after  spillage  of  crude 
oil  for  6  minutes.  On  March  10  the  fire  was  successfully  extinguished, 
and  oil  continued  to  spill  for  approximately  1  monLh  before  all  wells  were 
capped. 
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Our  observations,  except  for  those  kept  and  made  available  to  us  by 
Coast  Guard  personnel,  were  confined  to  two  periods,  March  5-11  and  March 
14-20. 

MEASUREMENT  TECHNIQUES 

Tides 

Tides  were  recorded  by  a  Hydro  Products  521  model  submerged  wave  and 
tide  gage.  This  hydrostatic  pressure-sensing  unit  is  designed  for  bottom 
installation,  but  since  there  was  little  time  to  dive  to  emplace  the 
instrument  properly,  it  was  simply  suspended  by  a  rope  from  the  MP-41M 
platform  in  firm  contact  with  the  bottom.  The  instrument  functioned  well 
and  gave  data  which  agree  very  well  with  tides  predicted  for  Pensacola 
(Fig.  4). 


Salinity-Temperature 

Salinity  and  temperature  were  determined  with  a  Beckman  RS-5  induction 


s alinometer.  The  induction  electrode  was  mounted  generally  just  above  the 


current  meter  and  raised  and  lowered  to  the  various  depths  at  which  currents 
were  monitored  so  that  simultaneous  temperature  and  salinity  readings  could 
be  made.  This  instrument  functioned  well  at  all  times.  A  few  values  above 


37  parts  per  thousand  are  anomalous  for  this  area  of  the  Gulf  of  Mexico 
and  probably  indicate  some  instability  in  calibration  or  misreading  of  the 
instrument .  Figure  5  is  a  record  of  salinity  and  temperature  for  the  entire 


period  of  observations  from  the  C.G.C.  Dependable . 

Winds 

Wind  speed  and  direction  were  recorded  at  quarter-hour  intervals  by 
the  men  on  watch  on  the  U.S. C.G.C.  Dependable  during  the  whole  period  of 
our  observations  on  board  that  vessel.  Crewmen  on  watch  on  the  cutters 
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Figure  4.  Observed  and  predicted  tides. 


Salinity  and  temperature  records  from  5-  and  10-foot  depths  at  Dependable  station 


Pt •  Spencer  and  Ft.  Estero  also  kept  logs  at  various  times  of  wind  speed 
and  direction  at  our  request.  The  lightship  New  Orleans  regularly  records 
meteorologic  data  every  6  hours  year-round  and  relays  this  information  *o 
the  Weather  Bureau.  Figure  2  gives  a  plot  of  wind  speed  and  one  of  wind 
direction  for  the  entire  month  of  March  based  on  the  lightship  records. 
Figure  6  shows  wind  speed  and  direction  recorded  from  March  15  to  March  20 
by  the  Dependable.  On  Figure  1  the  lightship  wind  records  for  March  are 
plotted  as  wind  rose  B  (based  on  115  observations)  for  direct  comparison 
with  the  winds  expected  in  an  "average"  March. 

Currents 

The  most  successful  current  monitoring  was  done  with  Marine  Advisers 
Q-15  ducted  current  meters,  which  were  raised  and  lowered  by  means  of  the 
ship's  winches.  Savonius  rotor-type  meters  were  employed  with  some  success 
at  fixed  levels  by  suspending  the  instruments  on  cables  from  the  cutter 
Reliance  and  from  the  lightship  New  Orleans . 

The  greatest  difficulty  in  measuring  currents  in  shallow  coastal 
waters  is  the  effect  of  the  oscillatory  orbital  wave  motions  on  the  speed 
sensors.  The  Savonius  rotor,  for  example,  will  spin  up  with  the  wave 
motions  and  can  overestimate  the  steady  current  by  a  factor  of  three  or 
four.  A  bidirectional  ducted  meter  such  as  the  Q-15,  used  extensively  in 
this  program,  largely  eliminates  this  error. 

The  chief  problem  in  current  monitoring  from  ships  is  the  motion  of 
the  ship  itself  even  when  at  anchor.  In  this  study  it  was  noted  that 
whenever  winds  exceeded  about  12  knots  the  vessels  began  to  yaw  from  side 
to  side,  swinging  on  the  anchor.  This  motion  of  the  ship  drags  the  current 


11 


Wind  tpnnd,  knots 


Morch,  197  0 


Figure  6.  Current  speed  and  direction  at  5-foot  depth  and  wind  speed 
and  direction  at  Dependable  station. 
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meter  and  causes  rhythmic  oscillations  of  readout,  giving  errors  that  can 
be  as  much  as  .5  knot  too  high.  This  type  of  error  is  obvious  on  the  record; 
all  erroneous  data  have  been  deleted  from  consideration.  Current  observa¬ 
tions  from  the  Dependable  are  plotted  on  Figures  6  and  8. 

The  vertical  motions  of  the  ship  are  an  additional  problem  which  can 
give  false  impressions  about  currents.  Both  side-to-side  swinging  and 
vertical  motion  led  to  a  generally  unsatisfactory  performance  for  the  Savonius 
meters,  and  information  from  these  is  not  used  here.  The  ducted  meters 
respond  in  a  selective  way  directionally.  For  a  detailed  discussion  of 
these  types  of  error  see  Murray  (1969) . 

Measuring  currents  very  near  to  the  surface  with  any  type  of  suspended 
meter  is  difficult  except  in  periods  of  calm  sea  states.  The  periods  of 
high  waves  encountered  in  the  study  sometimes  made  it  impossible  to  obtain 
current  readings  at  5  feet  depth  because  of  extreme  swinging  of  the  instru¬ 
ment  as  waves  rocked  the  ships.  Use  of  floating  drogues  or  drift  cards  for 
surface  current  measures  therefore  seems  preferable  if  adequate  means  of 
tracking  are  available. 

Sometimes  currents  changed  in  speed  very  rapidly,  as  can  be  seen  in 
Figure  6  at.  about  0700  CST,  when  current  speed  changed  from  .05  to  1.9  knots 
in  a  2-hour  period.  This  poses  a  sampling  problem.  The  intervals  for  com¬ 
pletion  of  a  single  profile  were  about  1  hour,  meaning  that  readings  at  any 
given  depth  were  repeated  about  once  an  hour.  With  such  rapid  changes  taking 
place,  this  may  not  have  been  sufficient  to  record  all  peak  speed  periods. 

It  would  be  preferable  to  have  currents  at  a  given  depth  monitored  con¬ 
tinuously  . 
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Ideally,  much  more  could  have  been  done,  and  will  be  done  as  the 
"stace  of  the  art"  for  shallow-water  current  measuring  Instruments  Improves. 
It  is  tempting  to  think  of  moored  strings  of  ducted  current  meters  spaced 
at  fixed  depths  ard  recording  continuously.  The  moored  deployment  would 
eliminate  many  of  the  problems  inherent  in  unsteady  ships,  the  continuous 
recording  would  give  better  resolution,  and  the  ducted  instruments  would 
give  better  directional  recponse.  Although  this  seeius  a  greatly  Improved 
means  of  profiling,  another  approach  will  be  necessary  for  continuous  study 
of  very-uear-surface  currents  at  times  of  varying  sea  state,  a  problem  of 
critical  Importance  in  the  movement  of  oil  on  water. 

Waves 

Wave  height  and  wave  period  were  recorded  by  the  Hydro  Products  521 
model  instrument  previously  described  in  the  section  on  tides.  As  mentioned, 
this  instrument  was  installed  by  suspending  it  from  a  platform.  Its  attach¬ 
ment  to  the  platform  probably  added  some  noise  to  the  data  (particularly  the 
wave  period)  owing  to  the  effect  of  the  structure  on  waves  approaching  from 
certain  directions.  Other  than  this,  the  records  seemed  excellent  and  are 
plotted  in  Figure  7  with  wind  speed  for  comparison. 

Choppy  seas  are  indicated  by  low  period  numbers.  At  1200  CST  on  March 
17,  for  instance,  during  a  high  wind,  wave  period  readied  a  minimum  of  3 
seconds.  Shortly  after  this  the  wind  fell  off  abruptly  and  the  record  began 
to  show  a  steady  swell  with  period  of  about  7  seconds.  The  swell  probably 
existed  prior  to  this  but  was  masked  by  the  seas. 

Slicks 

The  Coast  Guard's  routine  Situation  Reports  gave  direction  and  approxi¬ 
mate  width  and  length  of  the  oil  slicks  when  visibility  conditions  allowed 
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observation  from  helicopters  or  other  aircraft.  Aerial  remote  sensing  by 
photography  and  infrared  scanners  also  recorded  the  distribution  of  oil 
on  several  occasions.  The  infrared  (8-14p)  scanning  device  showed  the  form 
and  xtent  of  the  oil  slick  especially  well.  In  these  images  the  lighter 
films  of  oil,  which  formed  the  greater  part  of  the  area  of  the  slick, 
generally  appeared  darker  than  the  surrounding  water.  The  heavier  agglomer¬ 
ated  masses  of  oil,  which  occurred  in  patches  up  to  a  centimeter  in  thickness, 
showed  as  very  light  areas  on  the  scanner  imagery. 

ANALYSIS  AND  INTERPRETATION  OF  DATA 

During  the  12  days  at  the  site  of  the  Chevron  spill  our  research  group 
collected  a  tremendous  volume  of  data.  This  data  was  subsequently  augmented 
by  numerous  photographs  and  imagery  taken  by  several  government  agencies  and 
private  companies,  the  information  in  the  Situation  Reports,  and  assorted 
other  data.  Because  of  time  limitations,  in-depth  analysis  has  been 
restricted  to  the  characteristics  of  the  currents;  gross  relationships 
between  measured  current,  winds,  and  slick  movements;  and  theoretical 
slick  analysis.  We  have  arranged  most  of  ‘.he  observational  data  in  fig¬ 
ures,  and  this  section  will  be  mainly  a  discussion  of  the  salient  features 
exhibited  by  the  figures. 

Tides 

On  Figure  2  the  most  conspicuous  feature  of  the  tides  predicted  for 
March  is  the  biweekly  cycle  of  high-range  tropical  tides  and  low-range 
equatorial  tides,  which  are  characteristic  of  this  area  of  the  Gulf  of 
Mexico.  The  maximum  range  for  tropic  tides  is  about  2.5  feet,  and  the 
range  for  equatorial  tides  may  be  only  0.1  foot  or  less.  In  restricted 
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waters  strong  winds  can  increase  or  diminish  tidal  range  considerably.  The 
tides  are  largely  of  diurnal  character,  but  semi-diurnal  components  may  be 
evident,  particularly  during  equatorial  tides. 

Figure  4  shows  the  measured  tides,  with  predicted  curves  plotted  for 
comparison.  Agreement  between  the  two  is  excellent,  indicating  that  the 
tide  predictions  at  Pensacola  are  useful  for  the  area  of  study.  The  inter¬ 
val  of  March  6-11  spans  a  period  of  equatorial  tides  and  that  of  March  15-20 
covers  a  period  transitional  from  tropic  to  equatorial.  The  semi-diurnal 
components  are  evident  for  the  period  of  March  8-9  and  again  in  the  period 
March  17-20. 

Salinity  and  Tempt. rature 

Salinity  and  temperature  profiles  in  the  area  of  the  41C  platform  were 
quite  variable  in  the  upper  portions,  while  temperature  and  salinity  at 
depths  below  15  feet  were  rather  constant  at  about  16°C  and  35  parts  per 
thousand,  respectively .  At  times  the  temperature  profiles  were  virtually 
isothermal,  but  a  distinct  halocline  generally  existed  between  5  and  15 
feet  beneath  the  surface.  During  periods  of  high  winds  and  waves  the 
halocline  moved  downward,  reflecting  deeper  mixing  (see,  for  instance,  the 
salinity  records  from  the  Dependab le  for  March  17  cited  in  the  appendix). 

In  the  zone  above  the  halocline,  salinity  and  temperature  variations 
at  a  station  were  extreme  and  sudden.  Figure  9  shows  the  situation  which 
existed  at  the  Dependable  just  prior  to  and  during  the  first  NASA  over¬ 
flight.  A  mass  of  highly  turbid  water  moved  past  the  Dependable  at  approxi¬ 
mately  1330  CST,  and  temperatures  quickly  declined  by  2°  at  the  surface  while 
salinity  dropped  from  24  to  13  parts  per  thousand.  Figure  5  shows  the  record 
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of  this  event  as  well  as  that  at  the  10-foot  depth,  where  declines  of 
temperature  end  Balinity  are  not  as  great.  A  similar  passage  of  a  freshwater 
body  was  recorded  in  the  £t.  Spencer,  station  1,  observations  on  the  night 
of  March  5,  when  salinities  declined  abruptly  from  27  to  18  parts  per  thous¬ 
and,  then  slowly  returned  to  normal  levels  (see  appendix). 

These  low- salinity  surficlal  water  masses  usuallv  were  clearly  demar¬ 
cated  from  higher  salinity  waters  by  a  sharp  "tide  line"  or  interface. 

Such  interfaces  are  familiar  phenomena  in  the  Mississippi  Delta  area.  The 
interface  is  prominent  not  only  because  of  the  scrong  contrast  in  turbidity 
across  it  but  also  because  frequently  there  are  concentrations  of  floating 
debris  along  it.  Oil  likewise  tended  to  hang  up  along  it.  Slicks  of  great 
width  would  contract  to  narrow  bands  along  the  interface  and  abruptly  assume 
the  same  orientation  as  the  interface. 

The  fact  that  oil  tended  to  be  stranded  along  this  boundary,  together 
with  the  fact  that  only  a  small  amount  of  oil  reached  shore  from  this  large 
spill  and  this  only  on  Breton  and  Grand  Cosier  islands,  suggests  that  coastal 
lew-salinity  water  masses  may  have  acted  as  a  buffer  zone,  protecting  the 
marshlands.  River  water  discharged  from  the  numerous  passes  of  the  flooding 
Mississippi  must  also  create  a  hydraulic  gradient  leading  to  a  ret  offshore 
displacement  of  surface  waters. 

Winds 

Wind  rose  B  of  Figure  1,  based  on  115  observations,  shows  v.  us 
observed  by  the  lightship  New  Orleans  in  March  1970.  While  the  "average" 
March  depicted  In  A  shows  a  clear  predominance  of  winds  from  the  east, 
south,  and  southeast,  the  March  1970  record  differs  in  the  great  frequency 
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of  winds  from  the  northwest.  This  was  perhaps  a  fortunate  circumstance  in 
preventing  spread  of  oil  into  the  marshlands  to  the  west  and  northwest. 

The  wind  speeds  and  directions  indicated  in  Figure  2  show  the  erratic 
nature  of  winds  during  March.  The  two  periods  of  highest  winds,  on  March  8 
and  March  17,  occurred  during  the  periods  of  observations  and  added  variety 
as  well  as  difficulty  the  operations.  These  periods  are  discussed  in 
the  section  on  waves.  Average  wind  speed  for  the  whole  month  is  approxi¬ 
mately  12  knots. 

Figure  6  shows  winds  recorded  on  the  Dependable  for  March  15-20,  dur¬ 
ing  most  of  which  time  winds  from  the  southeast  prevailed. 

Waves 

Our  data  on  significant  wave  height  and  average  wave  period  near  41C, 
displayed  In  Figure  7,  are  restricted  to  the  periods  of  on-site  observa¬ 
tions,  March  6-11  and  March  15-20.  Significant  wave  height  is  the 

average  of  the  highest  one-third  of  the  waves.  The  correlative  wind  speed 
in  the  figure  shows  that  each  observation  period  had  a  clearly  defined  storm 
pattern.  During  the  lightest  winds  the  wave  height  diminished  to  0.2-0. 3 
foot,  while  the  highest  waves  ("4.0  feet),  early  on  March  18,  were  actually 
swell  generated  offshore.  The  expected  close  correlation  between  wind 
speed  and  wave  height  is  obvious,  but  a  close  examination  does  reveal  some 
interesting  features. 

During  the  first  storm  period  the  wind  blew  from  the  north,  a  direc¬ 
tion  in  which  the  fetch  is  limited  to  about  10  nautical  miles.  According 
to  Wiegel's  (1962)  empirical  equation,  the  maximum  possible  wave  height 
Hj^  given  in  a  nondimensional  form 
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where  U  Is  the  10-meter  wind  speed  and  F  is  the  fetch.  This  expression  pre¬ 
dicts  a  maximum  of  4  feet  for  but  the  wind  duration  was  not  sufficient 

to  generate  this  height  of  wave;  the  height  reached  only  2.7  feet.  In 
addition,  the  fetch  is  too  limited  for  any  appreciable  swell  to  develop, 
and  this  is  demonstrated  by  the  abrupt  decrease  in  wave  height  in  phase  with 
the  decay  of  the  wind  speed. 

The  second  storm  (March  16-18)  blew  from  the  southeast,  the  direction 
of  maximum  fetch  into  this  area.  The  wind  speed  attained  its  maximum  at 
1200  hours  on  March  17,  but  the  maximum  observed  wave  height  of  4.2  feet 
lagged  by  12  hours  because  of  the  arrival  of  large  swell  from  offshore. 

The  interaction  of  sea  and  swell  is  better  illustrated  by  the  wave 
period  diagrams  in  Figure  7.  During  the  first  storm,  as  soon  as  the  wind 
speed  exceeded  10  knots  tne  9ea  began  to  build  and  the  wave  period  dropped 
abruptly  to  5.3  seconds,  reflecting  the  surface  chop  (wind  waves),  then 
gradually  recovered  up  to  8  seconds  with  the  arrival  of  swell.  Again,  in 
the  second  storm,  the  sea  began  to  build  when  the  wind  exceeded  10  knots. 

The  wave  period  dropped  graduaJly  as  low  as  3  seconds  until  1200  hours 
March  17,  when  the  wind  speed  suddenly  dropped  from  30  knots  to  10  knots. 

This  discontinuity  is  immediately  reflected  in  the  wave  period,  with  a 
step-up  in  period  from  3  seconds  to  7  seconds  as  the  sea  quickly  dissipated 
and  the  swell  predominated. 

Prior  studies  (Johnson  and  Hwang,  1961;  Masch,  1961)  have  shown  that 
the  wave  height  and  period  can  significantly  affect  the  dispersal  of  fluid 
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properties  such  as  amount  of  suspended  solid  waste  or  thermal  maxima,  and 
so  by  extension  the  surface  oil  concentration.  It  is  unfortunate  that  with 
the  excellent  quality  of  the  wave  data  during  these  two  intervals  there  is 
practically  no  quantitative  data  on  the  oil  slick  with  which  to  make  com¬ 
parison.  The  slick  geometry  at  other  times  can,  however,  be  related  to  the 
wave  field  indirectly  through  the  wind  speed — this  in  discussed  in  the 
section  on  diffusion  theory. 

Currents 

A  priori  considerations  and  experience  lead  us  to  expect  a  mixture 
of  current  types  in  this  area— tidal  currents,  wind-drift  currents,  slope 
currents  from  setup,  setdown,  and  from  Mississippi  River  outflow,  and 
relative  currents  associated  with  freshwater-saltwater  density  gradients. 

Near-Surface  Currents. — The  unreduced  current  speeds  and  directions 
from  the  5-foot  depth  during  March  15-20  are  plotted  on  Figure  6,  together 
with  the  wind  speeds  and  directions.  (X/ing  mainly  to  exceptional  logis¬ 
tical  support,  our  most  reliable  current  information  comes  from  this  period. 
The  figure  shows  that  the  wind  remained  consistent  from  the  southeastern 
quadrant  but  that  the  current  direction  rotated  clockwise  in  tnree  broad 
cycles.  The  periods  were  usually  well  over  24  hours;  so  these  are  not 
simple  rotary  tidal  currents.  The  wind  speed  was  generally  greater  than 
10  knots  and  the  current  speed  generally  exceeded  0.6  knot. 

A  clear  relationship  between  the  win.,  and  the  5-foot  current  existed 
during  only  one  time  interval.  Between  1800  hours  March  16  and  1100  hours 
March  17  the  wind  blew  from  the  southeast  at  over  20  knots;  at  the  same 
time  the  current  flowed  into  the  northwesterly  quadrant  at  about  0.5  knot. 


The  Situation  Reports  disclose  that  the  only  severe  pollution  of  beaches  on 
the  adjacent  islands  took  place  on  the  afternoon  and  evening  of  March  16  on 
Breton  Island.  Oil  was  carried  onshore  by  the  currents  driven  by  the  strong 
southeasterly  winds. 

Since  the  magnitude  of  tidal  currents  is  approximately  proportional 
to  the  tidal  range,  the  diminishing  current  on  March  19  probably  reflects 
the  oncoming  equatorial  tide  of  March  20. 

The  overall  lack  of  correlation  between  the  wind  and  the  unreduced 
current  observed  at  a  depth  of  5  feet  suggests  that  the  other  current 
agents  operating  (tidal,  relative,  slope)  are  usually  of  the  same  magnitude 
or  greater  than  the  wind  current  contribution,  effectively  masking  any 
relationships  from  all  but  a  detailed  analysis.  It  seems  probable  now 
that  the  5-foot  current  does  not  reliably  represent  the  surface  current. 

It  is  unfortunate  that  the  vertical  ship  motion  and  wave  action  prevented 
us  from  getting  any  detailed  data  closer  to  the  surface  than  5  feet. 

Subsurface  Currents. — Two  further  levels  of  the  current  observations 
during  March  16-20  are  shown  in  Figure  8.  The  general  magnitude  of  the 
current  speed  differs  little  et  the  25-foot  level  and  the  45-foot  level 
from  that  of  the  near-surface  current.  Speeds  of  0.4  to  0.5  knot  are  most 
common  at  these  depths.  The  current  directions  at  the  25-  and  45-foot 
levels,  while  remarkably  similar  to  one  another,  differ  sharply  from  the 
near-surface  current  directions  in  that  they  are  reasonably  steady.  During 
most  of  this  observation  period  these  currents  were  moving  toward  the 
south-southwest.  There  is  a  possibility  that  these  seaward- flowing  cur¬ 
rents  were  returning  to  the  offshore  water  which  was  brought  onshore  by 


the  steady  southeasterly  wind.  This  seaward  return  flow  effect  at  depth  is 
quite  clear  on  March  16  and  17,  when  the  surface  waters  were  moving  onshore 
with  the  oil  that  washed  onto  Breton  Island. 

In  summary,  during  the  observation  interval  at  the  Dependable  station, 
changes  of  current  speed  and  direction  below  10  feet  are  minimal.  There  is 
also  ■  decided  reduction  in  the  tendency  toward  rotary  flow  in  the  lower 
layers. 

Measurements  (see  appendix)  from  the  C.G.C.  Pt_.  Spence r  in  40  feet  of 
water  on  March  5  and  6  were  made  during  the  calmest  period  of  the  month; 
they  clearly  showed  a  clockwise  rotating  diurnal  tidal  current  with  minimum 
interference  from  wind  currents.  The  current  flowed  initially  south  after 
high  tide,  turned  westerly,  and  then  flowed  north  about  13  hours  after  high 
tide.  The  average  speed  in  the  upper  10-foot  layer  was  1  knot;  below  this 
level  the  average  speed  was  0.5  knot.  Directional  turning  of  the  lower 
layer  lags  slightly  behind  that  of  the  upper  layer. 

Tidal  Currents, — In  order  to  isolate  end  examine  the  periodic  tidal 
currents  from  the  mixture  of  types  comprising  the  observed  current,  harmonic 
analysis  was  carried  out  on  the  5-  foot  level  currents  at  the  Dependable 
station  March  16-19  and  on  the  5-,  15-,  25-,  and  45-foot  levels  for  March 
18. 

Harmonic  analysis  extracts  the  periodic  components  which  can  be 
ac.cunted  for  by  any  set  number  of  harmonics.  The  first  two  harmonics, 
namely  the  diurnal  tidal  current,  with  a  period  of  1  lunar  day  (24  hours, 

50  minutes),  and  the  semi-diurnal  tidal  current,  with  a  period  of  12  hours, 
25  minutes,  account  for  nearly  95  percent  of  the  total  periodic  tidal  cur¬ 
rent.  As  shown  in  Figure  10 ,  we  nave  used  the  first  three  harmonics  in 
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our  analysis.  In  this  figure,  which  illustrates  the  method,  the  three 
diagrams  to  the  left  are  the  currents  accounted  for  by  the  three  harmonics 
individually — the  irregular  diagram  on  the  right  is  the  resultant  current 
obtained  by  adding  all  the  harmonics  together.  Such  diagrams  are  called 
tidal  ellipses;  the  current  vector  at  any  hour  is  given  by  an  arrow  drawn 
from  the  origin  to  the  point  labeled  with  that  hour  number.  The  length 
and  direction  of  the  arrow  give  the  current  speed  and  direction,  respec¬ 
tively.  On  the  graphs  only  the  current  vector  at  hour  1  is  drawn.  The 
smooth  curve  only  serves  as  an  envelope  to  connect  the  end  points  of  the 
arrows.  The  semi-diurnal  tide  goes  around  its  ellipse  twice  each  day,  etc. 
As  the  diurnal  and  semi-diurnal  tide  waves  usually  differ  in  magnitude, 
direction,  and  phase,  the  resultant  tidal  ellipse  may  be  a  rather  compli¬ 
cated  figure.  As  expected  from  our  knowledge  of  the  tides  in  the  Gulf  of 
Mexico,  the  diurnal  tidal  current  is  the  most  dominant.  Because  of  the 
brevity  of  these  records,  i t  is  not  possible  to  assign  an  absolute  time 
to  the  hourly  notations  on  the  tidal  ellipses. 

With  the  effects  of  steady  currents  removed,  the  tidal  ellipses  of 
Figure  11  show  the  day-to-day  change  in  the  tidal  currents  at  the  b-foot 
level.  The  direction  of  rotation,  as  it  should  be  in  the  northern  hemi¬ 
sphere,  is  clockwise  owing  to  Coriolis  forces.  There  is  an  obvious  decrease 
in  the  magnitude  of  the  tidal  current  as  time  progresses  because  of  the 
fact  that  we  are  nearly  at  equatorial  tide  on  March  19  (the  strength  of 
the  current  is  1.0  knot  on  March  16  and  only  0.6  knot  on  March  19).  The 
orientation  of  the  ellipse  is  generally  in  an  onshore-offshore  direction — 
the  maximum  currents,  then,  are  onshore  and  offshore  with  the  flood  and 
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ebb  of  the  tide.  This  analysis  has  successfully  removed  a  fairly  typical 
tidal  current  from  the  rather  confused  data  presented  in  Figure  6. 

The  changes  in  the  tidal  current  ellipse  c.s  a  xunction  of  depth  are 
shown  in  Figure  12.  Increasing  frictional  effects  as  the  bottom  is 
approached  considerably  distort  the  tidal  current  ellipse,  resulting  in 
occasional  periods  of  counterclockwi se  rotation.  Speeds  are  lower  near 
the  bottom  and  the  higher  speeds  have  a  shorter  duration. 

It  is  apparent  that  oil  purposely  sunk  into  the.  bottom  layers  during 
equatorial  tides  would  have  the  least  aid  from  the  tidal  currents  in  dis¬ 
persing. 

Slicks 

Figure  13,  which  is  a  plot  of  wind  vectors  against  the  slick  vectors 
observed  by  helicopter  visual  surveillance,  gives  the  distinct  visual 
impression  that  there  is  a  degree  of  correlation  of  wind  direction  and 
slick  direction.  A  more  detailed  examination  of  the  angle  between  simul¬ 
taneous  wind  and  current  vectors  is  not  possible  owing  to  (1)  uncertainties 
in  the  time  of  slick  observations  and  (2)  the  fact  that  the  only  detailed 
wind  data  are  restricted  to  the  Dependable  observations  of  March  15-20. 

It  does  seem  clear  from  the  figure,  however,  that  after  the  wind  has 
blown,  say,  a  day  from  a  given  direction  the  slick  aligns  with  the  wind. 

This  relationship  implies  a  good  correlation  between  the  movement 
of  the  upper  foot  of  surface  water  and  the  wind.  We  have  already  shown 
that  there  is  no  clear  relation  between  the  current  at  the  5-foot  level 
and  the  wind.  Figure  14  shows  the  velocity  profile  predicted  by  Reid 
(1959)  for  a  pure  wind-drift  current.  There  is  a  thin  high-velocity 
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'  Figure  14.  Non-dimensional  vertical  velocity  frofile  for  pure  wind- 

drift  current  by  Reid  (1957).  Z/D  is  relative  iepth,  U  as  current 
speed,  t  is  fluid  density,  and  t,  is  shear  stress  uniform  throughout 
the  fluid. 
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layer  predicted  near  the  surface;  it  is  this  type  of  situation  which  would 
allow  the  very-near-surface  water  to  move  with  the  wind  while  the  under¬ 
lying  waters  (say,  at  5-foot  depth)  are  only  barely  affected.  Some  other 
means  of  monitoring  currents  closer  to  the  surface  is  necessary  to  clarify 
the  problem. 

At  times  of  high  winds  with  steady  direction,  oil  did  definitely  move 
with  the  wind,  as  was  discussed  in  the  section  on  currents  with  respect  to 
the  strong  winds  of  March  16-17,  which  pushed  oil  onto  Breton  Island.  The 
U.S.  Geological  Survey  photomosaics  and  color  photographs  of  oil  slicks 
which  were  provided  to  us  showed  five  cases  in  which  oil  slicks  extended 
with  the  wind  (within  20  degrees),  as  indicated  by  plumes  of  spray  from 
the  fire  boats  hosing  down  the  platform  (Fig.  15).  Two  cases  showed  no 
clear  indication  of  wind  direction,  and  one  showed  a  slick  extending  at 
about  90  degrees  to  the  wind  direction. 

As  mentioned  previously  .  L’tk  orientation  often  changed  abruptly 
at  interfacial  lines.  Slicks  ;>lsw  changed  strikingly  in  shape  at  inter¬ 
facial  lines,  becoming  narrow  bands  (Fig,  16)  of  thicker  oil  accumulation, 
with  coalescing  "rafts"  or  "mats"  of  oil  up  to  1  centimeter  thick.  This 
change  in  shape  of  the  slick  and  the  observed  stranding  of  other  floating 
objects  along  the  interface  suggest  that  the  boundary  is  a  locus  of  down- 
welling  at  which  floating  materials  collect.  The  slicks  seem  literally 
to  contract  to  lie  along  the  interface  as  if  "swept"  into  a  linear  accu¬ 
mulation.  The  narrow  accumulations  are  also  sometimes  distorted  in 
patterns  suggestive  of  very  large  Kelvin-Helmholtz  waves,  which  may  form 
along  the  interface  of  fluids  in  motion  relative  to  each  other. 
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Figun  X  .1  row  band  of  oil  stranded  along  brackish 

water-saltwater  interface.  (Provided  by  U.S.G.S. ) 
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Figure  15.  Photograph  oi  slick,  March  14,  1970.  (Provided  by  U.S.G.S 


Barge  line  was  for  trapping  onshore-flowing  oil.  For  direction'  and 


scale  refer  to  Fig.  IS.) 
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Wltui'e  inter!  actis  were  not  Involved,  slicks  showed  shapes  that  were 
straight  or  arcuate  to  varying  degrees.  The  arcuate  forms  are  indicative 
ot  currents,  either  tidal  or  wind  generated,  which  rotated  as  the  slick 
developed.  Straight  forms  develop  during  periods  of  steady  current  direc- 
t  ion. 

The  "mat"  type  of  oil  accumulation  and  the  thin  film  slicks  were  the 
only  types  of  floating  oil  noted.  Some  reports  were  heard  at  the  time  of 
thicker  accumulations,  but  they  were  not  substantiated  by  our  observations. 

The  most  useful  images  of  oil  slicks  that  were  available  for  the 
study  were  those  made  by  8-14u  infrared  scanning  imagery.  These  showed 
the  thin  film  as  dark  areas  relative  to  the  surrounding  waters  and  the 
"rafts"  of  oil  in  thicker  accumulations  as  bright  areas.  The  same  imagery 
showed  excellent  delineation  of  interfaces  between  different  water  masses 
an  well.  Sensitivity  of  this  type  of  imagery  at  high  altitude  allowed 
sensing  of  the  entire  slick  on  a  single  image  with  adequate  detail  of 
size  and  shape.  Future  emergencies  of  this  kind  should  be  better  docu¬ 
mented  by  this  type  of  sensing. 

THEORETICAL  CONSIDERATIONS— SLICK  PREDICTION 

Introduction 

The  mixing,  dispersion,  or  diffusion  of  oil  on  the  sea  is  undoubtedly 
a  complicated  process  under  field  conditions  because  the  process  is  largely 
controlled  by  changes  in  the  weather  and  sea  states.  However,  our  under¬ 
standing  of  the  physical  processes  and  ultimately  of  the  spread  of  the  oil 
can  be  greatly  increased  by  application  of  the  theory  of  turbulent  diffusion. 
This  approach  has  been  successfully  applied  in  the  diffusion  of  atmospheric 
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pollutants  and  in  studies  of  the  spread  of  dye  and  other  substances  on 
lakes  and  on  the  sea. 

Simple  turbulent  diffusion  theory  does  not  consider  such  effects  as 
particle  fallout,  particle  interaction,  coagulation,  or  evaporation.  Spe¬ 
cial  terms  may  be  added  to  the  equations  to  account  for  these  effects, 
however  (see  Yudine,  1959;  McNown  and  Lin,  1952;  and  Sutton,  1955).  Empiri¬ 
cal  determination  of  the  diffusion  coefficient  K  in  the  field  during  an  oil 
spill  should  take  these  effects  into  account  to  some  degree.  The  following 
development  applies  only  to  an  oil  slick  dispersing  in  a  laterally  homo¬ 
geneous  water  body.  As  shown  in  the  previous  section,  horizontal  density 
gradients  and  discontinuities  will  considerably  distort  the  diffusion 
pattern. 

Statistical  Theory  of  Turbulence 

In  his  pioneering  paper  of  1921  and  again  in  a  detailed  treatise  in 
1935,  G.  I.  Taylor  developed  the  statistical  theory  of  turbulence  which 
laid  the  groundwork  for  most  modern  studies  of  diffusion  problems  in  the 
atmosphere  and  ocean.  By  considering  the  one-dimensional  spread  of  par¬ 
ticles  (or  marked  fluid  parcq)>.  from  a  continuous  point  source,  Taylor 
developed  the  relation  for  the  y  direction 

_  T  t 

Ry  (O  d£di  (1) 

o  o 

2  2 
where  o  is  the  variance  of  the  fluid  parcels  around  the  source,  v’  is 

y 

the  square  of  the  turbulence  intensity,  assumed  to  be  constant  over  a  period 
of  a*  *  ral  hours,  v'  -  v  -  v  as  usual,  £  and  t  are  time  parameters,  and 
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R(0  is  the  Lagrangian  autocorrelation  function.  For  a  detailed  derivation 
of  (1)  see  Haltiner  and  Martin  (1959,  p.  156).  The  Lagrangian  autocorrela¬ 
tion  function  measures  the  relationship  between  the  velocity  of  a  parcel  at 
one  instant  and  che  velocity  of  the  same  parcel  at  some  subsequent  instant. 


(2) 


Hence  for  very  short  time  lags  v'(t)  =  v'(t~0  and  R^(0  -  1.  For  long  time 
lags  there  will  be  on  the  average  no  relation  between  v'(t)  and  v'(t~0;  so 
R^(0  -  0.  Similar  expressions  for  and  R^(0  extend  the  Taylor  analysis 

to  three  dimensions. 


After  considering  very  short  times,  performing  the  double  integration, 
and  then  taking  tl-.e  first  derivative  with  respect  to  time  of  ,  equation 
(1)  reduces  to 


-  (v,2)1/2  .  (3) 

V  * 

In  the  case  of  long  diffusion  times  /  R  (£)  dE,  *  t  a  constant,  and  equation 

j0  y  y 

(1)  reduces  to 

.2  * 

-  2  v'*  t  .  (4) 

The  subscript  *  rc  ars  to  the  ultimate  constant  value  this  quantity  will 
attain.  It  a  uniform  horizontal  mean  flow  is  passing  the  source  moving  in 
the  positive  x  direction  at  a  speed  U  assumed  constant  over  a  period  of 
several  hours,  the  relation  T  ■  x/U  can  bf.  used  to  transform  the  spread  of 
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the  material  into  the  x-y  plane;  so  (3)  becomes  (for  short  diffusion  times) 


do 

dx 


(5) 


and  (4)  becomes  (for  long  diffusion  times) 


do 


dx 


2  v' 2  £ 


(6) 


*  * 

where  i  the  Lagrangian  eddy  size  ■  U  t 


Thus  4^-  for  short  times  and 
dx 


do 

dx 


for  long  times  are  both  constants.  Figure  17  illustrates  the  three 


ensuing  diffusion  regimes:  (a)  the  linear  or  short-term  regime,  where 

1/2 

a  «  x,  (b)  the  long-term  regime,  where  o  “  x  ,  and  (c)  the  intermediate 

y  y 

or  transitional  regime,  where  “  xn,  1/2  <  n  <  1.  For  the  purposes  of  the 
illustration  the  are?  occupied  by  the  dispersing  material  can  be  outlined  by 
assuming  that  the  radius  to  the  visible  edge  of  the  slick  (y  *  R)  can  be 
identified  initially  with  some  constant  multiple  of  the  root-mean-square 
displacement  c^.  It  will  also  be  assumed  that  the  visible  radius  corres¬ 
ponds  to  an  equal  concentration  C  contour.  Pritchard  (1966)  has  expressed 
objection  to  this  assumption  in  dye  studies  of  diffusion  processes,  but 
because  of  the  technical  difficulty  of  obtaining  cross  slick  measurements 
of  the  concentration  C  the  assumption  is  at  present  unavoidable.  An  expand¬ 
ing  slick  should,  then,  initially  exhibit  a  linear  growth  of  its  radius, 
followed  ultimately  by  a  slow  parabolic  expansion.  We  can  logically  apply 
the  one-dimensional  Taylor  equations  to  an  expanding  oil  slick  because 
(a)  oil  is  generally  constrained  by  buoyancy  forces  to  the  surface  (no 
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Linear,  transitional,  and  parabolic  diffusion  regimes  predicted  by  Taylor  theory 


mixing  in  the  vertical  direction)  and  (b)  the  longitudinal  diffusion  of  the 
oil  particles  in  the  direction  of  the  mean  horizontal  velocity  is  negligible 
compared  to  the  displacements  in  this  direction  by  the  mean  current  itself. 
Taylor  has  further  demonstrated  that  the  diffusion  coefficient  K  can 

y 

be  expressed  as 


Survey  and  Remote  Sensing  Inc.  of  Houston,  clearly  and  repeatedly  shows  the 
presence  of  the  linear  regime  near  the  spilling  well,  41C.  This  is  excep¬ 
tionally  well  shown  in  Figure  18.  If  we  assume  first  that  c  =  ■<.  R  (the 
root-mean-square  displacement  is  directly  proportional  to  the  visible  radius 


R)  and  as  a  first  approximation  that  the  proportionality  constant  a  is  near 
unity  (0^  :  R) ,  we  can  use  the  photographs  and  imagery  to  obtain  estimates 
of  the  relative  turbulence  intensity  of  the  oil  from  (5)  and  the  diffusion 


coefficient  of  the  oil  from  (7). 


Table  1  lists  the  source  and  values  for  the  relative  turbulence 
2.  1/2 

intensity  (v'  )  /U  of  the  oil,  measured  as  the  tangent  of  one-half  the 
angle  of  expansion  20  of  the  oil  slick  from  the  well. 
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Table  1 


Intensity  of  Turbulence  Indicated  by  Oil 
Source _ tan  £  =  (v'2  )l!2 /M 


U.S.G.S 

.  Infrared  Color,  undated 

0.41 

U.S.G.S 

-  Color 

00074,  3/14 

0.25 

U.S.G.S 

.  Color 

00003,  3/16 

0.  33 

U.S.G.S 

.  Color 

00071,  3/14 

0.28 

U.S.G.S 

.  B  6  W 

Mosaic,  3/22 

0.22 

U.S.G.S 

.  B  &  W 

Mosaic,  3/13 

0.22 

Remote  1 

Sensing 

Inc. ,  3/14 

0.22 

Remote  Sensing 

Inc. ,  3/11 

0.10 

Average 

0.25 

These  results  are  encouraging,  as  it  is  well  known  (see,  e.g.,  Kalinske  and 
Pien,  1944;  and  Bowden,  1962)  that  the  relative  turbulence  intensity  in 
channel  flow  is  in  the  range  of  0.05  to  0.15.  The  effect  of  wave  motions 
on  the  oil,  especially  breaking  action,  could  certainly  increase  this  value 
up  to  the  observed  average  0.25.  Although  the  data  points  are  few,  a  plot 
of  tan  0  against  wind  speed  in  Figure  2  does  seem  to  indicate  a  higher 
relative  turbulence  intensity  of  oil  at  higher  wind  speeds  (stronger  wave 
action) . 

Unfortunately,  there  is  only  one  set  of  imagery  which  shows  a  slick 
with  sufficient  length  and  contrast  to  permit  the  measurements  of  the 
parabolic  regime  and  so  the  determination  of  a  value  for  the  diffusion 
coefficient  K^.  This  is  the  Remote  Sensing  Inc.  infrared  imagery  ( 8 — 1 4 u 
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scanner)  ot  March  II,  1230  hours  CST.  The  plot  of  against  x  is  shown 

in  Figure  19.  As  predicted  by  (6),  the  linear  relationship  is  strong;  the 

2 

regression  coefficient  “  0.97.  With  d'  /dx]^  now  known,  use  of  (7)  leads 

5  2 

to  an  estimate  of  the  diffusion  coefficient  of  oil  X  of  4  x  10  cm  /sec  at 
that  time.  Since  Figure  7  allows  that  wind  speed  and  wave  height  were 
moderate  on  March  11,  4  x  10^  cm^/see  is  probably  a  low  rather  than  a  high 
estimate.  Comparison  with  the  few  published  measurements  of  horizontal  eddy 
viscosity  (the  diffusion  coefficient  for  momentum)  is  encouraging.  Sverdrup, 
Johnson,  and  Fleming  (1942)  report  values  (p.  485)  ranging  from  2  x  10^ 
cm  /sec  in  a  small  area  with  weak  currents  to  4  x  10  cm  /sec  in  strong 
currents.  It  is  now  generally  understood  that  the  diffusion  coefficient 
is  a  function  of  the  scale  of  the  phenomena  under  investigation.  As  Orlob 

(1959)  clearly  points  out,  the  boundary  of  the  field  (basin),  among  other 

* 

parameters,  will  restrict  the  maximum  eddy  size  that  may  develop;  and 

*  2  * 

since  is  related  to  by  two  other  constants  *  (v1  ) /U  (8), 

the  diffusion  coefficient  in  a  bounded  area  should  likewise  approach  a 
constant  value. 

It  should  be  pointed  out  that  the  y  coordinate  of  as  one  moves  in 
the  positive  direction  (downslick)  does  not  necessarily  follow  an  equal 
concentration  contour.  The  evidence  in  Figure  19  does,  however,  seem  to 
Indicate  that  this  is  approximately  the  case  for  some  distance  downslick, 
until  the  concentration  level  at  a  distance  o  from  the  average  slick  center 

y 

line  C(y  -  a  )  falls  below  the  concentration  level  which  is  assumed  to  serve 


as  a  visible  limit.  In  Figure  19  the  discontinuity  in  the  oy  versus  x 
plot  at  x  :  30,000  feet  is  Interpreted  as  that  point  where  first 
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Figure  19.  Plot  of  slick  width  variance  (o  )  against  distance  from 
source  (x)  ,  illustrating  linear  relation  that  determines  the  diffusion 
coefficient. 


45 


C(y  -  j  )  <  C (y  -  R). 

Despite  Its  power  and  apparent  applicability  to  oil  slicks,  the  Taylor 
theory  cannot  predict  concentration  contours  or  slick  boundaries.  For  this 
we  must  turn  to  the  classical  diffusion  equations. 

Genjeral  Diffusion  Equations  and  Slick  Prediction 

Again  regarding  the  spread  of  oil  as  a  two-dimensional  process  in  the 
x-y  plane,  with  a  near  horizontal  motion  U  in  the  postive  x  direction,  the 
diffusion  equation  is 


3C 

3t 


<Kx  + 


f-  (K 

3y  y 


ac. 

9yJ 


If  the  slick  has  come  to  steady  state  —  *  0,  and  assuming  (1)  that  the 

O  t 

diffusion  in  bitropic  K  -  K  **  K  and  (2)  that  K  has  reached  its  constant 
r  x  y 

value,  (9)  becomes 


(9) 


,g. ,(*.&)■ 

For  a  derivation  of  the  complete  diffusion  equation  see  Haltiner  and  Martin 
(1959,  p.  265).  The  exact  solution  of  (10)  for  a  two-dimensional  point 
source  (infinite  line  source)  in  an  infinite  fluid  is 


C 


Q  eU  x/2K 
2ttK 


K 

o 


(ID 


where  Q  is  the  rate  of  emission  of  the  source  (g/sec),  Kq(A)  is,  using 

standard  notation,  the  modified  Bessel  function  of  2nd  Kind  and  Zero  Order, 

-»  2  1/2 

and  r  *  (x“  +  y  )  .  For  very  reasonable  values  of  U  and  K  the  first  term 
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of  the  asymptotic  expansion  for  Kq(‘)  (Abramowitz  and  Stegun,  196b,  eq. 
9.7.2;  Roberts,  1923,  p.  644)  can  be  substituted  in  (11),  yielding 


Ux/2K. 


TUeret ore 


Ur 

q—  (1  -  COS  i) 


2  (?K  U  r) 


where  if  is  measured  in  a  clockwise  direction  from  the  positive  x  axis.  Equa¬ 
tion  (13),  then,  is  the  expression  we  are  seeking;  it  gives  the  distribution 
of  the  emitted  oil  as  a  function  of  distance  and  direction  (r,  1)  from  the 
source,  the  rate  of  emission  at  the  source  Q,  the  ambient  current  speed  U , 
and  the  oil  diffusion  coefficient  K.  The  value  of  Q  (personal  communication. 
Commander  Dixon,  U.S.C.G.)  during  the  Chevron  incident  is  taken  as  1000  bbl/day 
1241  g/sec. 

Furthermore,  following  Roberts  (1923)  it  can  easily  be  shown  that  the 
maximum  half-width  of  a  given  concentration  contour  is 


ymax  .  .1/2 

(2-re)  LC 

and  by  setting  $  »  0  in  (33)  the  maximum  extension  downslick  of  a  given  con¬ 
centration  contour  is 


max  ,  „,J 
«»T,  KuC 


It  is  not  difficult  conceptually  to  relate  the  oil  concentration  C  in 


g/cm  to  the  slick  thickness.  There  may,  in  fact,  be  a  direct  relationship 


if  the  hydrostatic  equation  p  *  agh  is  used.  If  C  is  identified  with  p  (as 
they  have  the  same  dimensions),  then  the  slick  thickness  h  is  directly 
proportional  to  C. 


This  idea,  however,  needs  much  closer  examination. 

One  difficulty  in  using  (13)  to  predict  slick  boundaries  is  that  it 
predicts  small  but  finite  oil  concentrations  at  infinitely  great  distances 
from  the  emitting  source.  In  meteorology  is  routine  (Sutton,  1953)  to 

take  10  percent  of  the  maximum  concentration  as  a  visible  boundary  contour. 
Unfortunately,  C  -*■  *  at  the  origin  in  (13);  so  we  must  take  some  concentra¬ 
tion  close  to  the  origin  as  a  control  level.  Inspection  of  the  solution  of 
(13)  along  the  center  line  of  the  slick  ( <f>  m  0)  for  K  =  4  x  10^  cm^/sec,  Q 
1241  g/sec,  and  U  «  50  cm/sec  (1  knot)  shows  that  the  presence  of  the  source 
at  r  =  0  is  just  beginning  to  affect  the  solution  adversely  in  the  vicinity 
of  r  *  100  meters.  The  control  concentration  for  the  Remote  Sensing  Inc. 

(RSI)  slick  of  March  11,  with  K  »  4  x  10"*  cm^/sec,  Q  *=  1241  g/sec,  and 

_3  2 

U  “  50  cm/sec,  Is,  using  (13),  0.7828  x  10  g/c.m  .  The  slick  boundary 

-  4  2 

C  n,  then,  is  taken  as  0.7828  x  10  g/cm  .  If  this  latter  value  for 
v  «  K 

the  concentration  at  the  boundary  is  uacd,  (14)  predicts  one-half  the  slick 
width  and  (15)  predicts  the  slick  length.  The  predicted  maximum  wid^h  of 
the  RSI  March  11  slick  is  1,534  meters,  the  observed  maximum  width  is 
1,507  meters;  the  corresponding  predicted  length  is  9,995  meters,  while 
the  observed  length  is  14,710  meters.  The  slick  width  is  predicted  more 
accurately  than  the  length,  but  the  results  indicate  the  basic  soundness 
of  our  approach.  Furthermore,  the  average  of  23  observations  on  slick 


L 


(16) 
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length  made  by  Coast  Guard  helicopter  visual  surveillance  teams,  reported 
in  the  Situation  Reports  and  displayed  in  Figure  13,  comes  out  to  12,600 
meters,  compared  to  our  prediction  of  9,995  meters. 

An  exact  expression  for  the  area  of  a  slick  is  not  easily  obtained 
from  the  basic  equations;  so  we  have  resorted  to  an  elliptical  approximation 


A  “  nab 


(17) 


where  a  the  semi-major  axis  is  one-half  the  s.'ick  length  given  by  (15)  and 
b  the  semi-minor  axis  is  the  half  width  given  by  (14).  Therefore 


A 


8  /2re  K  V 


R 


(18) 


Inspection  of  (14)  shows  that  the  slick  width  will  increase  directly 

with  Q  and  decrease  directly  as  U  increases.  It  is  noteworthy  that  (14) 

predicts  that  the  slick  width  is  independent  of  K.  The  slick  length  is 

more  strongly  dependent  on  Q,  increasing  as  and  decreasing  directly  as 

3 

U  or  K.  increases.  The  area  given  by  (18)  increases  as  Q  and  decreases  as 
2 

U  and  the  first  power  of  K. 

The  effect  of  varying  U,  with  Q  and  K  held  constant,  is  shown  in 
Figure  20.  Low  values  of  U  produce  a  long,  large  slick;  values  of  U  approxi¬ 
mately  to  those  observed  in  oie  Held  during  th-  Chevron  maueia 

(0.6-0. 8  knot)  produce  relatively  short  and  slender  slicks  much  like  those 
repeatedly  observed  in  photographs  and  visually.  The  estimated  time  t 
required  to  reach  a  steady  state  is  also  shown  in  Figure  20.  This  is  crudely 
approximated  as  the  time  required  for  oil  to  travel  with  the  current  out  to 
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Figure  20.  Comparison  of  observed  slick  against  theoretical 
current  velocity  on  slick  size. 


the  outermost  slick  boundary.  A  numerical  solution  of  the  time  dependent 
differential  equation  should  produce  a  considerably  11101“  accurate  estimate. 

The  shape  factor  SF  of  the  slick,  defined  as  the  ratio  of  slick  width 
to  slick  length,  is,  combining  (14)  and  (15), 


SF 


4  (2 *)1/2  K  C 


Q  « 


1/2 


R 


It  is  noteworthy  that  the  relative  shape  of  che  slick  is  independent  of  the 
current  speed,  being  controlled  solely  by  (K/Q)  for  a  given  boundary  con¬ 
centration.  Low  values  of  (K/Q)  produce  elongate  slicks;  high  values  of 
(K/Q)  produce  more  circular  slicks.  For  the  Chevron  incident,  using 

c  2  __  )  2 

K  »  4  x  10  cm  / oec,  n  R  =  .78  x  10  g/cm  ,  and  Q  =*  1241  g/sec,  the 

predicted  value  of  the  shape  factor  SF  *■  0.15.  Table  2  lists  the  shape 
factors  determined  from  slick  geometries  reconstructed  from  U.S.  Coast 
Guard  helicopter  visual  surveillance  Situation  Reports.  Unfortunately,  only 
data  from  the  last  week  of  the  incident  were  of  sufficient  detail  to  be  of 
use  in  this  regard. 

Table  2 

Shape  Factors  During  Chevron  Incident 


Date 

Time 

Shape  Fac 

3/26 

0011  Z 

0.50 

3/26 

1756  Z 

0.19 

3/29 

0140  Z 

0.22 

3/29 

1905  Z 

0.07 

3/30 

1521  Z 

0.17 

3/30 

1940  Z 

0.16 

3/31 

203C  Z 

0.32 

Average 

0.23 
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T!n:  appro*  lmate  agreement  of  the  average  SF  with  the  predicted  value  of  0.15 
suggests  that  the  value  of  K  on  the  average  should  be  closer  to  6  x  10^  than 
to  4  x  1Q'>  cnT/sec  for  the  last  week  of  March. 

Figures  21  through  27  present  in  more  detail  and  in  graphic  form  most 
of  the  results  we  have  discussed  up  to  now.  Figures  21,  22,  23,  and  24  show 
slick  area  as  a  function  of  Q,  K,  and  U,  whereas  Figures  25,  26,  and  27  show 
slick  width  and  slick  length  as  functions  of  these  same  three  variables.  In 

further  studies  along  these  lines  an  effort  should  be  made  <.o  incorporate 
these  parameters  into  a  single  nondimens ional  diagram  or  nomogram. 

Although  this  two-dimensional  diffusion  approach  seems  to  present  a 
good  approximation  of  gross  slick  geometries,  it  is  really  based  on  the 
concentration  value  at  the  slick  boundary,  for  which  we  made  only  a  rough 
estimate  (10  percent  of  C  at  r  -  100  m,  $  *  0).  In  future  work  this  point 
must  be  examined  in  much  greater  detail. 

SUMMARY 

Observations  prior  to  and  during  the  long  period  of  continuous  oil 
spillage  from  a  Chevron  production  platform  near  the  Mississippi  Delta 
permit  certain  generalizations  to  be  made  about  movement  of  oil  in  coastal 
shallow  water  bodies.  An  understanding  of  the  situation  seemed  to  require 
information  on  currents,  winds,  tides,  waves,  and  orientations  and  shapes 
of  spreading  oi.1  masses.  Salinity  ana  temperature  structure  of  the  sur¬ 
rounding  waters  seemed  useful  for  indicating  the  riverine  and  marine 
influences  expected  in  the  delta.  Information  was  gathered  on  all  of  these 
in  the  hope  that  understanding  could  be  gained  of  the  principal  agents  which 
control  rates  and  direction  of  motion  of  oil  on  water. 
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Figure  24.  Slick  area  (A)  as  a  function  of  diffusion  coefficient  (K) 
for  various  values  of  current  speed  (U) . 
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Figure  27.  Slick  width  (W)  and  slick  length  (L)  as  a  function  of 
diffusion  coefficient  (K)  for  various  values  of  current  speed  (U) 


Winds  during  March  were  somewhat  more  varied  than  would  be  expected 
from  long-term  records.  Strong  winds  to  the  southeast,  which  are  usually 
rather  infrequent  in  March,  were  probably  important  in  reducing  the  amount 
of  oil  reaching  shore. 

At  lease  once  each  week  winds  were  sufficiently  strong  to  generate 
waves  of  2  to  4  feet,  significant  wave  height.  Two  such  periods  of  wave 
buildup  and  decay  were  documented  by  our  instrumentation.  There  is  no 
doubt  that  the  stirring  and  mixing  caused  by  these  waves  played  a  sig¬ 
nificant  role  in  the  continued  breakup  and  diffusion  of  the  oil  slicks. 

Salinity  and  temperature  in  the  upper  15  feet  of  the  water  column 
showed  rapid  and  strong  variation.  These  were  related  to  masses  of  turbid 
lower  salinity  water  deriving  from  the  river  mouths.  Interfaces  between 
the  river-derived  and  marine  water  masses  were  usually  sharp  and  acted 
as  barriers  to  movement  of  oil  by  other  agents.  The  large  amount  of  water 
being  discharged  from  the  Mississippi  during  March  was  one  of  the  more 
important  deterrents  to  movement  of  oil  into  the  marshes  because  of  this 
barrier  effect,  plus  the  escape  of  discharged  river  water  out  to  sea. 

Observations  on  the  local  tides  indicated  that  they  agreed  very  well 
with  those  predicted  at  Pensacola.  The  current  observations,  on  the  other 
hand,  showed  a  complicated  pattern;  although  tidal  components  were  strong, 
their  interaction  with  wind-drift  currents  and  slope  currents  masked  any 
clear  relationship  '*■  en  the  tide  and  the  current  or  the  wind  and  the 
current.  Harmonic  na3ysis  of  the  currents  did  isolate  typical  current 
ellipses  which  h«d  strong  diurnal  and  semi-diurnal  components  (0. 5-1.0  knot). 
This  analysis  further  showed  tnat  the  amplitude  of  the  tidal  current  decreased 
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with  decreasing  tidal  range  (as  the  time  of  equatorial  tides  approached) 
and  increased  with  distance  from  the  bottom. 

A  general  association  of  slick  movement  with  wind  direction  is  indi¬ 
cated  from  aerial  photographs  and  from  Coest  Guard  Situation  Reports.  Oil 
washed  ashore  on  Breton  Island  on  March  16  and  17  in  a  period  of  strong 
winds  blowing  from  the  well  to  the  island.  This  was  the  rain  period  when 
oil  reached  shore.  The  indicated  association  of  slick  movement  with  wind 
is  in  contrast  with  the  lack  of  correlation  between  the  wind  and  the  cur¬ 
rent  at  the  5-foot  level.  One  possible  explanation  suggested  by  tneory 
Is  that  the  slick  was  embedded  in  a  thin,  high-velocity  surface  driven 
by  the  wind  stress. 

Slick  shapes  and  orientations  were  often  controlled  by  long  (up  to 
several  miles)  lines  or  zones  separating  brackish  water  masses  from  those 
--of  normal  Gulf  salinity.  There  apparently  was  no  horizontal  component  of 
motion  across  these  interfaces,  as  oil  and  other  material  frequently 
accumulated  along  them.  These  discontinuities  often  formed  a  natural 
barrier,  preventing  oil  from  moving  onto  the  shore. 

A  theoretical  analysis  was  undertaken  to  ascertain  the  effects  of 
the  rate  of  oil  discharge,  the  ambient  current  speed,  and  the  rate  of  oil 
diffusion  or.  the  area,  width,  and  length  of  an  ideal  slick.  Typical  values 
for  the  Chevron  incident  predict  an  equilibrium  or  steady-state  slick 
length  of  about  10  kilometers  in  length  and  1.5  kilometers  in  width. 

These  values  agree  very  well  with  the  on-site  observations.  Other  results 
on  the  theoretical  analysis  are  as  follows:  (1)  For  a  given  oil  spill 
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incident  the  current  speed  largely  determines  the  size  a  slick  will  reach; 
(2)  the  ratio  of  the  diffusion  coefficient  to  the  rate  of  oil  discharge 
controls  the  shape  of  the  slick;  (3)  high  discharges  will  provide  elongate 
slicks;  and  (4)  the  consistently  high  current  speeds  during  the  Chevron 
incident  kept  the  slick  size  to  a  minimum  and  greatly  enhanced  the  rate 
of  diffusion  of  the  oil. 

CONCLUSIONS 

1.  In  the  waters  surrounding  the  Chevron  MP-41C  platform  it  is 
apparent  that  wind  stress  tends  to  drive  a  thin  layer  of  surface  water 
before  the  wind.  This  is  indicated  by  the  frequent  agreement  of  slick 
orientation  with  wind  direction. 

2.  Detailed  measurements  of  currents  at  5-foot  depth  did  not  show 
agreement  with  slick  orientations  or  wind  directions  except  for  periods  of 
higher  wind  speeds. 

3.  Harmonic  analysis  removed  significant  diurnal  and  semi-diurnal 
components  from  the  current  records.  The  remaining  non-periodic  components 
largely  lect  a  combination  of  wind-generated  currents  and  slope  currents 
from  setup,  setdown,  and  river  discharge. 

4.  Tidal  currents  showed  a  strong  dependence  on  the  range  of  tide, 
with  the  speed  dropping  by  a  factor  of  2  or  3  from  tropic  tide  periods  to 
equatorial  tide  periods.  Dispersion  of  oil  could  be  expected  to  diminish 
accordingly. 

5.  Currents  at  depth  showed  a  greater  importance  of  non-periodic 
components;  and,  although  rotational  effec  .  were  still  evident,  the  deeper 
currents  frequently  ran  counter  to  the  near-surface  currents.  The  tra¬ 
jectory  of  chemically  dispersed  oil  could  therefore  not  be  predicted  from 
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surface  observations. 

6.  The  sharp  boundaries  associated  with  lower  salinity  water  masses 
moving  through  the  area  were  sites  of  oil  stranding  and  accumulation. 

7.  Theoretical  analysis  indicates  a  strong  effect  of  current  speed 
on  slick  size.  The  high  speeds  of  the  Chevron  incident  rapidly  diffused 

the  oil  and  kept  slick  sizes  to  a  minimum.  Slick  size  and  shape  also  varied 

as  a  function  of  rate  of  oil  discharged  and  the  coefficient  of  diffusion, 

5  2 

which  was  empirically  determined  to  be  4  x  10  cm  /sec  for  the  Chevron  case. 


ACKNOWLEDGMENTS 


Our  gratitude  is  expressed  to  the  Geography  Programs,  Office  of  Naval 
Research,  and  to  the  U.S.  Coast  Guard  for  enabling  us  to  carry  out  this 
study.  In  particular,  cur  thanks  go  to  Miss  Evelyn  L.  Pruitt,  ONR;  Admiral 
Ross  T.  Bullard,  Eighth  Coast  Guard  District  Commander;  Lieutenant  Commander 
David  H.  Dickson,  Jr.,  Chief,  Intelligence  and  Law  Enforcement  Branch  (Eighth 
District);  Captain  Lewis  W.  Tibbies,  Jr.,  Chief  of  Staff,  Eighth  U.S,  Coast 
Guard  District;  Lieutenant  Commander  Roy  L.  Foote,  Operations,  Eighth  Dis¬ 
trict;  and  Lieutenant  Commander  Stephen  L.  Richmond  of  Operations,  New 
Orleans  Group. 

We  are  grateful  for  the  assistance  of  the  officers  and  crews  of  the 
U.S.C.G.  cucters  Dependable ,  Reliance ,  Pt.  Spencer,  and  Pt_.  Estero.  Special 
thanks  go  to  BMCM  Rufus  R.  McGiness  and  BMCM  Clyde  M.  Bowden. 

Photographic  coverage  of  the  oil  slicks  was  provided  by  the  U.S. 
Geological  Survey  through  the  assistance  of  Messrs.  William  A.  Radlir.ski 
and  Charles  J.  Rovinove.  Mr.  Jerry  Wermund  of  Remote  Sensing  Inc. ,  Houston, 
Texas,  provided  excellent  infrared  scanner  imagery  of  the  slicks. 

The  capability  of  the  Coastal  Studies  Institute,  which  permitted  a 
rapid  response  to  this  emergency  situation,  was  created  by  the  long-term 
support  of  the  Institute  by  the  Office  of  Naval  Research,  Geography  Programs, 
under  Contract  No.  N00014-69-A-0211-0003 ,  Project  No.  NR  388  002. 

Mr.  Norwood  Rector,  Mr.  Stan  Sandifer,  Miss  Sharon  Dunn,  and  Mr. 

William  E.  Hart  of  the  Coastal  Studies  Institute  contributed  long  hours 
beyond  the  call  of  duty  to  assure  successful  completion  of  the  project. 


64 


APPENDIX 


r 


&r! 

u> 

in 

o 

in 

O 

O 

m 

<n 

!X> 

CM 

sf 

m 

m 

m 

m 

cn 

'O 

■‘T 

m 

00 

00 

m 

r*N 

m 

m 

cn 

cn 

r~> 

rH 

M 

G; 

in 

<N 

o 

i.n 

CM 

in 

Mi¬ 

m 

CM 

rH 

O 

m 

o 

m 

cn 

CM 

rH 

rH 

O 

CM 

rH 

o 

m 

<r 

<}• 

CM 

CM 

a 

-» - 

r< 

CM 

cm 

CM 

rN 

r— ■J 

-H 

m 

en 

cn 

m 

rn 

cn 

m 

CM 

rH 

o 

O 

O 

O 

g 

o 

S 

tM 

CM 

CM 

rH 

rH 

rH 

r-H 

f-«  ■ 

cm 

CM 

CM 

CM 

CM 

<M 

CN 

cm 

CM 

CM 

<M 

CM 

CM 

CM 

(M 

o 

o 

o 

o 

o 

o 

O 

O 

O 

o 

o 

o 

O 

o 

1 

\ 


si 

•rH. 

CM 

<N 

O 

O 

n 

O 

o 

MJ 

lO 

00 

in 

M? 

rH 

O 

i£> 

00 

00 

o 

\0 

O 

vO 

*•* 

vO 

00 

CM 

i 

O I 

M> 

in 

M3 

m 

lO 

lO 

M3 

CO 

CO 

rH 

H 

O' 

m 

CM 

O 

r-*. 

00 

rH 

O' 

O' 

00 

r-H 

CM 

<x> 

r-H 

cn 

O' 

MD 

■ 

^ ; 

r — ^ 

rH 

cm 

rn 

cn 

m 

rH 

« — c 

CM 

CM 

rH 

CM 

CM 

cn 

i — i 

rH 

rH 

CM 

rH 

rH 

rH 

CM 

CM 

rH 

CM 

CM 

CM 

rH 

*H 

sj 

■! 

Si 

<S) 


u 


>, 


v> 


O  CT  in  O  O  O 

O  r>  in  h  h  h  fsj 


m  m  o  in  in  o  •/> 

(NrNiO^HnfNH 


U">  iT\  iT>  O  O 
r'-*  m  <n  CM 


00 

•H 


m  mm  o  c0  m  o 

N  fN  h>  IT)  cn  n  t”!  fS 


n  n  m  in  n  q>  m 


m 

•vj\OsjninciNC0 


OOtXJCOvOOCOCOl/l 

COHHhniOOOu) 


mommoofNON 

COOOMmOO\£MO 


vo  io  vo  \o  in  m 


\0\0\0^^miinin 


m  .-h  r>»  o  O'  rv  ■srr^min<T<NO^ 


mjvOvOvO'CiAlAiA 


oimocoNHmH 


iA\jD\oso\c»ninin 

hhhhhhhh 


o  O'  O'  <r  o  o  <r 
noNHnnmio 


H  H  fNJ  ><r  sj  >■}  lT| 

m  m  cn  cn  m  m  cn 


NOfN'jsrminmi 
ch  cn  m  n  m  n  n  tn 


cN^-cM-s-Hrmmm 

NMmmnmnn 


co  o  O'!  <r  m  m  m 

rH  n  m  rn  m  n  cn 


in 


m 


& 


cm  m  O  m  O  m  m 

H  H  N  CM  m 


cm  m  o  m  o 

r-H  r-H  CM 


O  O 
cm  cn 


m 

m 


m 

cNmomomo'O 

h  h  fM  cn  m  m 


in 

cMinomomoin 

H  H  <M  CM  ^  m 


O 

u 

Cm 


CM 


cn 


>9* 


C 

o 


CTJ 

4J 

V) 


u 

(0 

a  h 


C 

<u  c 

Cu  O 

CO  -H 


u  u 

On  CO 


S 

T 


65 


i 


Station  Date _ Profile  Depth  Salinity  Temperature _ Current  Speed _ Direction _ Tips 

5  2.5  21-46  16.15  1.5  261  0345 
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Oceanographic  observations  near  the  Chevron  spilling  well  off  the  Mississippi  Delta 
in  March  1970  revealed  relative  roles  of  various  physical  factors  of  the  regional 
estuarine  system  in  the  behavior  of  oil  slacks .( ^Surface  stress  from  the  wind  was 
most  important,  at  speeds  above  15  mph  the  slick  orientation  was  generally  determined 
by  the  wind  direction.  The  wind  also  indiractly  Effected  oil  which  was  sunk  by  dis¬ 
persant  in  the  wind  waves  promoted  mixing,  which  in  turn  affected  the  vertical  stabi¬ 
lity,  hence  eventually  the  velocity  profile.  Wind  setups  and  setdowns  were  correlated 
with  downward  and  upward  isopycnal  movements,  respectively.  Both  calculations  and 
observations  showed  that  tidal  currents  produced  an  L-shaped  slick  geometry  when 
winds  were  below  about  15  mph.  The  diurnal  rotation  of  the  tidal  currents  served  to 
limit  the  excursion  length  of  oil  from  the  source,  keeping  it  short  of  the  nearest 
shore.  The  presence  of  fresh  water  from  the  Mississippi  River  in  the  surface  layer 
and  the  consequent  development  of  convergence  lines  often  formed  a  natural  barrier, 
preventing  oil  from  encroaching  upon  the  shore.  Theoretical  analysis  using  turbulent" 
diffusion  theory  disclosed  that  the  area  and  length  of  a  steady-state  oil  slick  in¬ 
creased  with  oil  discharge  but  decreased  with  current  speed  and  the  lateral  diffusion 
coefficient.  Slick  width  increased  with  the  rate  of  discharge  and  decreased  with 
current  speed  but  was  independent  of  the  diffusion  coefficient.  The  width/length 
ratio  of  an  oil  slick,  which  was  notably  independent  of  current  speed,  was  controlled 
only  by  the  diffusion  coefficient  and  the  oil  discharge.  In  order  to  be  able  to 
predict  oil  slick  behavior,  understanding  of  the  interaction  between  these  multiple 
factors  is  essential  and,  for  that  end,  the  use  of  a  multiple-sensor  array  with  auto- 

D0  »'Onu  1473  matic  data  transfer  and  processing  capability  is  recommended. 

/  Security  CUs»ific«:ion 


44444 


